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ABSTRACT 

We study a class of early dark energy (EDE) models, in which, unlike in standard dark energy models, a 
substantial amount of dark energy exists in the matter-dominated era. We self-consistently include dark energy 
perturbations, and show that these models may be successfully constrained using future observations of galaxy 
clusters, in particular the redshift abundance, and the Sunyaev-Zel'dovich (SZ) power spectrum. We make 
predictions for EDE models, as well as ACDM for incoming X-ray (eROSITA) and microwave (South Pole 
Telescope) observations. We show that galaxy clusters' mass function and the SZ power spectrum will put 
strong constraints both on the equation of state of dark energy today and the redshift at which EDE transits to 
present-day ACDM like behavior for these models, thus providing complementary information to the geometric 
probes of dark energy. Not including perturbations in EDE models leads to those models being practically 
indistinguishable from ACDM. An MCMC analysis of future galaxy cluster surveys provides constraints for 
EDE parameters that are competitive with and complementary to background expansion observations such as 
supernovae. 
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1. INTRODUCTION 

Over the last decade, observational evidence has mounted 
in favor of dark energy, the mysterious component which 
dominates the energy content of the universe at present, 
and causes the expa n sion of the univers e to accelerate 
dKowalski et al. 20081: iKomatsu etal. 2009b . The nature 
of dark energy is one of the most tantalizing mysteries 
of present day cosmology. The simplest model for dark 
energy, the cosmol ogical constant, fits the current data well 
(Kowalsk i etal. 20081) . however, there are no strong con- 
straints on the time evolution of dark energy at present. Thus, 
evolving models of dark energy remain as alternative candi- 
dates for dark energy. Many non-cosmological constant mod- 
els have been suggested for dark energy, including scalar field 
quintessence models, modificati ons of the Einstein frame- 
work of grav it y, etc. (see rev iews Sahni & Starobinskv 2000: 
Carroll 200 it [Peebles & Ratra 20031: ICopeland et al. 20061 
Noiiri & Odinstov 20071: " ISahni & Starobinskv 2 006: 
Frieman et al. 20081 and references therein), but as of 
now there is no clear consensus on the nature on dark 
energy. An interesting class of models which have been 
suggested in the literature are early dark energy (EDE) 
models, in which the early universe contained a substantial 
amount of dark energy. Geometric probes of dark energy put 
strong constraints on the present-day nature of dark energy, 
and these constraints are expected to improve with future 
surveys. However, very little is known as to the nature of 
dark energy at early times due to paucity of data beyond 
redshifts of few. It has thus been hypothesized that even 



if dark energy at present behaves like the cosmological 
constant, in the past it could have had completely different 
behavior, leading to the idea of the EDE models. Dif- 
ferent fa cets of these models ha ve been studied in recent 
works dDodelson et al. 20001 : ISkordis & Albrecht 20021: 
iDoran & Robbers 2006b and references therein, and have 
been analyzed with respect to observa t ions extensively in 
recent times in ([Ti nder & Robbers 2008; Fran cis et al. 2008: 
Grossi & Springel 2009: Fedeli et al. 2009; Xi a & Viel 2009: 



Jennings et al. 2010; Alam 2010; de Putter et al. 2010). Cur- 



rent data places some constraints on these models, but does 
not rule them out. In this work we use a parameterization of 
the equation of state of dark energy to study the possibility of 
constraining EDE models using future large scale structure 
surveys. 

Ground and space-based telescopes targeting clusters 
in microwave [A tacama Cosmology Telescope (ACT) 
dHincks et al. 2009b. the South Pole Telescope (SPT) 
( Stan iszewski et al. 2009b and Planck l ] have begun opera- 
tion, or will return data shortly. Also current (Ebeling et 
al. 2010) and future missions have been planned for detect- 
ing clusters in X-ray waveband (Predehl et al. 2007). It has 
been recognized for some time now, that galaxy cluster sur- 
veys in X-ray or microwave via Sunyaev-Zel'dovich (SZ) ef- 
fect could be precision probes of cosmological parameters, 
in particul ar the dark energy de nsity and its equation of state 
parameter (lHaiman et al. 2001b . This is because redshift dis- 
tribution of massive clusters is exponentially sensitive to the 
growth of structure history of the universe, which in turn, 
bears the signature of dark energy. The presence of early dark 
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energy reduce structure formation, consequently the cluster 
abundances reduce. Thus cluster counts can provide a smok- 
ing gun probe for detecting early dark energy. For example, 
(Man tz et al. 20 10) have used a flux limited X-ray cluster data 
to constrain early dark energy component assuming the red- 
shift of transition between to 1 . The current generation of 
SZ surveys will also measure the power spectrum of the cos- 
mic microwave background with an unprecedented accuracy 
down to scales of 1 arc-min. The SZ surveys have started mea- 
suring the SZ power spect rum which is the do minant signal at 
scales of few arc-min ([H lFowler et al. 2010|) . The amplitude 
of the SZ power spectrum is proportional to the total number 
of objects that have formed and hence is extremely sensitive 
to the amount of dark energy present at early epoch. 

The paper is organized as follows. In Section [2] we explain 
the early dark energy formalism and model used for this anal- 
ysis. Section [3] expounds the nature of X-ray and microwave 
observations and mass-observable scaling relations. In sec- 
tion [4] we show predictions from current and future observa- 
tions, comparing our EDE models to a fiducial ACDM cos- 
mology. Finally, section [5] is devoted to conclusions and dis- 
cussion. 

2. EARLY DARK ENERGY 

Dark energy perturbations for dynamic dark en- 
ergy models have been studied in a number of works, 
usually under the formalism of a minimally coupled 
scalar fi e ld (See Maet al. 19991: iHwang & Noh 2001 



Hu 20021: iMakmarti & Liddle 20021: IWeller & Lewis 2 003. 
Bean & Pore 20041: iDutta & Maor 20071: iMota etai lOOl: 
Novosvadlvi & Sersiienko 2008; Jassal 20091 and refer- 
ences therein). In th is work we follow the formalism of 
( Well er & Lewis 2003[) . First order perturbations in a homo- 
geneous and isotropic large scale universe described by the 
Friedman-Lemaitre-Robertson- Walker (FLRW) metric take 
the form 

ds 2 = a 2 (r]) [(l+2^(x,r ? ))^ 2 -(l+2<l>(x,7 ? ))(5 a/ 3(ix a ^] , 

(1) 

where r\ is the conformal time, x is the length element, a(r\) is 
the scale factor, and <£, ^ are the Bardeen potentials. If proper 
isotropy of the medium is zero, then = 

Along with the matter and radiation components, we con- 
sider dark energy to be an additional fluid component, so that 
the dark energy perturbations are characterized by an equation 
of state and an adiabatic sound speed- 
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Defining the frame-invariant quantity c 2 si (the fluid sound 
speed in the frame comoving with the fluid), the evolution 
equations for a fluid with equation of state w,- = Ptj Pu and 
adiabatic speed of sound c 2 a i = pij p t , can be written as (prime 
denotes derivative with respect to rj) 



S[ = -3H(c 2 1 - wdSi - 9H 2 (c 2 1 - c 2 a ,)( 1 + w t ) y 



-(1 + wdkvi - 3(1 +WiW 



v' i =-U(l-3c 2 Jv i + 



kc lA 
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(4) 
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where T~L = a' /a is the conformal Hubble parameter and A is 
the acceleration (A = in the synchronous gauge, A = -\P in 



the Newtonian gauge). Adiabatic initial conditions are con- 
sidered. For the matter component, w m = c 2 = c 2 = 0. For the 
dark energy component, a fluid with varying wde > — 1 has 
c 2 de = wde- [dwoE I d(ln a)] / 3(1 +wde)- For scalar field like 
dark energy models, c 2 DE = 1. For a more general class of 
models, such as k-essence, c 2 DE could be variable as well. 

If we consider dark energy without perturbations, the quan- 
tities 5de , S'de - 0' anc * me matter density contrast is given by 

S f ^-HS f m - 47rGp m S m = . (6) 

thus the dark energy component appears only in the damp- 
ing term, so that a non-negligible amount of dark energy 
would lead simply to a suppression of clustering of matter 
at large scales. Not taking into account the dark energy per- 
turbations can le ad to gauge dependent results, as shown in 
( Park et al. 20091) . in this paper we have used the commonly 
used comoving gauge. If dark energy perturbations are in- 
cluded self-consistently, the results are gauge-independent. 

To study EDE models under this formulation, we consider a 
w-parameterization whi ch may represent a large class of vary- 
ing dark energy models ( Corasaniti et al. 20 03) 



w(a) = wo + (w m - w ) 



l _ e -(a-D/A t 
\+ e -(a-a t )/A t l_ e l/A t 



(7) 



where wo is the equation of state of dark energy today, w m 
is the equation of state in the matter dominated era, a t is the 
scale factor at which the transition between wq and w m takes 
place, and A t is the width of the transition. For studying EDE 
models with this parameterization, we choose w m > -0.1, to 
ensure the presence of adequate amount of dark energy at 
early times. 

Current data puts im pressive c onstraints on the values of 
these parameters (see Ala m 2010|) . however, many interest- 
ing EDE models still fall in the acceptable range. We choose 
two such EDE models in order to study the possibility of con- 
straining these models further using future observations. The 
models chosen are 

• Model 1 : w =-1.0, w m = -0.05, a t = 0.17, A t = 0.17; 

• Model 2 : w =-0.9, w m = -0.05, a t = 0.1, A, = 0.1. 

Model 1 behaves like ACDM at present, while Model 2 
has a higher value of the equation of state at present. 
We compare the results for these models with those for a 
ACDM model with identical values for the non-dark en- 
ergy parameters. The non-dark energy parameters are cho- 
sen from the WMAP7 CMB+BAO+HST best-fit, e.g., Vt h h 2 = 
0.0226, Q c h 2 = 0.1123,/* = 0.704, ^ = 0.963, as well as other 
parameters such as t he scalar amplitude A s and the reioniza- 
tion optical depth r dKomatsu et al. 20091). We modified the 
publicly available CQSMQMC code dLewis & Bridle 20021) 
for generating the transfer functions, and therefore the power 
spectrum at different redshifts for these models. 

Figure Q] shows the behavior of the equation of state of dark 
energy, and the matter power spectrum at redshift z = for 
these two models. For comparison we also plot the results for 
the matter power spectrum when dark energy perturbations 
are not considered. We see that when the perturbations are not 
considered, the EDE matter power spectrum, although there is 
a suppression the matter power spectrum, this effect is negli- 
gible, and for both models the matter power spectrum is very 
close to that for ACDM. However, when the perturbations are 
considered, the matter power spectrum is significantly differ- 
ent from ACDM, resulting in a ag = 0.74 for EDE1, and a 
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Fig. 1 . — Equation of state (panel (a)) and matter power spectrum at z = (panel (b)) for EDE models with wq = -1 .0, w m = -0.05, a t = 0. 17, A f = 
0.17 (EDE1, red lines), and with wo = -0.9, w m = -0.05, a t = 0.1, A t - = 0.1 (EDE2, blue lines). The black line in each panel represents the 
corresponding ACDM model, in panel (b), the dashed lines represent EDE without DE perturbations, while the solid lines represent EDE with 
perturbations. 

cr 8 = 0.76 for EDE2, whereas the ACDM has a <r 8 = 0.81 . 
This behavior is typical of EDE models with rapid and large 
transition of the equation of state. This translates to a rapid 
change in the dark energy perturbations at high scales (low k), 
thus to enhancement in the transfer function at low k. When 
normalized to the CMB scale, this effect shows up as an ap- 
parent suppression of power at high k, and therefore a low 
dg. The effect is stronger in EDE1 (i.e. lower cr 8 ) because the 
difference between the equation of state at present and in the 
matter dominated era in this case is larger. Further details on 
how the change in equation of state affects the d ark energy 
perturbations can be found in paper I (lAlam2010h . 

3. OBSERVATIONS 

Paper I (|Alam 2 010) explores effects of early dark energy 
on CMB power spectrum and other current observations; here 
we focus on forthcoming large-scale structure probes. We 
consider the influence of EDE on cluster counts and SZ power 
spectrum. 

3.1. Cluster Counts 

Redshift evolution of cluster abundance provides a valu- 
able insight into global dynamics of the universe. This abun- 
dance is sensitive to both expansion and growth history of the 
universe, and as such can be used to constrain standard cos- 
mological parameters like ag or ^om (see, e.g. Haiman et 
al. 2001), but also to explore possible mechanisms for ob- 
served acceleration, like modifying the law of gravity. 

For a spatially flat cosmologies considered here (Q K =0), 
all-sky number of objects more massive than M min in a redshift 
bin is: 



-10 
-12 




LCDM 
EDE 1 
EDE 2 



10 ! 



1 2 
Redshift 



1 2 
Redshift 



N(z) = 



47rr 2 (z)dr(z) dn(M,z) 



dz 



dM 



dzdM , (8) 



where the first term is the comoving volume, and the second 
is the differential mass function. 

From numerical simulations, we have a good quantification 
of the halo mass function. Although its precision is still not 



Fig. 2. — eROSITA survey, from upper left panel in clockwise direc- 
tion: minimal detectable mass (dashed line shows instrument capability, solid 
line shows minimal mass for galaxy clusters only with mass threshold of 
1O 14 M0), luminosity distance, comoving volume element covered by the sur- 
vey, and comoving cluster abundance as a function of redshift. 

at the level needed to make it insignificant for future cos- 
mological parameter studies (Cunha & Evrard 2009; Wu et 
al. 2009), that will not be relevant here as we want to explore 
the relative difference between our EDE models and the fidu- 
cial ACDM cosmology. The important point, on which we 
rely here, is that the mass function can be presented in a uni- 
versal form, where the cosmology dependence comes from 
the linear power spectrum, and the linear growth factor. At 
a 20% level (often better than that), it is indeed proven to be 
the case, for many different cosmological models (Jenkins et 
al. 2001; Linder & Jenkins 2003; Jennings et al. 2010; Bhat- 
tacharya et al. 2010) and for a wide range of redshifts (Heit- 
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mann et al. 2006; Lukic et al. 2007; Reed et al. 2007). Since 
we are interested in modeling X-ray and SZ observations, we 
will use mass function of halos defined as spheres enclosing a 
given overdensity, defined with respect to the critical density 
for the closure, p c (for analysis of different halo definitions 
see, e.g. White 2002; Lukic et al. 2009). 

One way of detecting groups and clusters of galaxies is 
via weak lensing maps (Marian & Bernstein 2006); this ap- 
proach is appealing as masses are probed directly there, but 
on a downside the method can probe only the most mas- 
sive systems, and only within limited redshifts (as galaxy 
shapes are very hard to measure beyond z ~ 1). Promising 
alternatives to the weak lensing are mass measures through 
Sunyaev-Zel'dovich effect (e.g. Carlstrom et al. 2002), X- 
ray emission (e.g. Kravtsov et al. 2006), and galaxy richness 
(High et al. 2010). For these, indirect mass measurements, 
one first has to connect cluster mass to a relevant obser v- 
able for a given survey, like X-ray flux (Mantzetal. 201 
or integrated Com pton y-parameter (IVanderlinde et al. 201 
Sehg al et al 2010|) . In this paper we will focus on X-ray and 
SZ surveys, and in the following we explain scaling relations 
we use. 

In the self- similar theory (Kaiser 1986, 1991), temperature 
of the intracluster medium scales with gravitational poten- 
tial (T oc M/R), and enclosed mass is M ex R 3 . This self- 
similar scaling was confirmed in hydrodynamical simulations 
(Mathiesen & Evrard 2001; Borgani et al. 2004; Kravtsov et 
al. 2006), as well as in observations (Vikhlinin et al. 2009). In 
the following, we will use M-T relation from Mathiesen & 
Evrard (2001): 



M .-l= 1Ol5 (?^) E{Z) ~ 



M 2 oo 



E(z)=H(z)/H . 

(9) 

We use the above relation, as it relates M 2 oo to the bolometric 
luminosity. While X-ray scaling relations are much tighter 
for M500 (Kravtsov et al. 2006), here we want to use the same 
mass for both X-ray and SZ estimates, and whether it is the 
best observational approach is not of concern here. 

Self-similar theory fails to correctly predict T — L relation 
(Markevitch 1998; Allen & Fabian 1998), and thus mass- 
luminosity relation is inaccurate as well. The departure from 
self- similarity is due to excess entropy in cluster cores which 
prevents gas from being compressed to very high densities 
(Ponman et al. 1999, Finoguenov et al. 2002). Therefore, in 
the following we will use M-L relation from Bartelmann & 
White (2003): 



L= 1.097 x 10 45 



200 



10 15 M^/r 



rE(z) 



1.554 



erg/s 



(10) 



X-ray observations are sensitive in a given energy band, 
whose flux is: 

Aird L (zr 

where is the band correction and diii) is luminosity dis- 
tance to a cluster: 



d L (z) = 



Ho Jo 



dz' 
E(z f ) 



(12) 



Emission from the intracluster medium is predominantly ther- 
mal bremsstrahlung, but for T < 2keV line emission from 
metals (clusters have metallicity ~ O.3Z ) becomes non- 
negligible (Sarazin 1986). Thus, to calculate we use 
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Fig. 3. — South Pole SZ survey, from upper left panel in clockwise di- 
rection: minimal detectable mass (dashed line shows instrument capability - 
same as in Figure |2j, angular diameter distance, comoving volume element 
covered by the survey, and comoving cluster abundance as a function of red- 
shift. 

XSPEC package (Arnaud 1996) and we assume Raymond- 
Smith (1977) plasma model for cluster emission. As the band 
correction depends on the plasma temperature, we solve it- 
eratively for minimum mass and corresponding temperature 
using XSPEC and equation^ Galactic absorption is modeled 
with constant column density of nu - 10 21 cm~ 2 , roughly cor- 
responding to Galactic mid-latitudes. In observational analy- 
sis one would use appropriate nu for the line of sight to each 
cluster, and column density will generally be smaller than our 
number for high latitudes, and larger for low latitudes. 

For the Sunyaev-Zel'dovich observations, one has to relate 
Compton y parameter (integrated over the solid angle) to the 
mass of an object. Furthermore, it is necessary to renormalize 
such relation to the directly observable SZ flux. In Fedeli et 
al. (2009), this was done using Sehgal et al. (2007) ^200-^200 
relation, and assumption that SZ signal outside cluster virial 
radius can safely be neglected, resulting in relation: 



6.763 x 10 



7 mJy / 



M 



200 



(d A (z)/lMpc) 2 V10 15 M G 



1.876 



\f{v)\E 2 '\z) . (13) 



Here, d A is angular diameter distance to the object, and f(y) 
is spectral signature of the thermal SZ effect: 

jcV r , 

x 4 



(e x -l) 2 



e x -\ 



hv 



(14) 



and T 7 is the CMB temperature. 



3.2. Sunyaev-Zeldovich Power Spectrum 

The current generation of CMB experiments will measure 
the power spectrum of the cosmic microwave background 
with an unprecedented accuracy down to scales of 1 arc-min. 
While the primary CMB fluctuations dominate the power 
spectrum at a degree scale, at scales of few arcminutes the 
secondary fluctuations arising from the SZ effect and lensing 
of the CMB becomes the dominant signal. 

Predictions for the SZ power spectrum amplitude C^sz 
(henceforth Q) can be made using the halo model and es- 
timates of the radial pressure profile of intra-cluster gas. 
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Assuming that the cluster gas resides in hydrostatic equi- 
librium in the potential well of the host dark matter halo, 
(Ko matsu & Se ljak 2002) demonstrated that the ensemble av- 
eraged power spectrum amplitude Q has an extremely sensi- 
tive dependence on ag where Q ex ag(^/z 2 ) 2 . The kinetic SZ 
power spectrum contribute roughly 10% to the total SZ angu- 
lar power spectrum and hence we consider only the thermal 
component of the SZ power spectrum (tSZ hereafter) here. 
The cosmological information of the tSZ power spectrum 
comes from the total number of halos of mass M > 1O 13 M / h 
in the survey area. Thus the SZ power spectrum contains a 
significant contribution from low mass clusters and group- 
mass objects. Note that measurements of the mass func- 
tion usually probe only massive halos (M > 1O 14 M ) and not 
group size halos. 

The SZ power spectrum can be calculated by simply sum- 
ming up the squared Fourier-space SZ profiles, y(M,z,f) of 
all clusters: 

<w./*£/*.^W (i5, 

where V(z) is the comoving volume per steradian and n(M, z) 
is the number density of objects of mass M at redshift z. g v 
is the frequency factor of the SZ effect. In this study, we 
show the power spectrum at 150 GHz where g y = -1. Note 
that whilst this calculation assu mes that clusters are Pois- 
son distributed, dShaw et al.2009|) have shown that including 
halo clustering modifies the power spectrum (compared to the 
Poisson case) by less than 1%. 

The number density of halos n(M, z) can be calculated for 
a ACDM cosmolog y using the fitting fun ctions provided by 
(|Jenkins et al. 200 lb and more recently (Tin ker et al. 2 008). 
For the SZ profiles (also mass and redshift dependent), pre- 
vious studies have frequentl y used either a simple /3 -model 
or the hydrostatic model of dKomatsu & Seljak 20021) . Note 
that both models simply assume that the gas resides in hy- 
drostatic equilibrium in the potential well of an NFW-like 
halo (and is isothermal in the case of the /3-model. Nei- 
ther model accounts for the fraction of hot gas that will have 
cooled and been converted into stars or a ny non-thermal en- 
ergy input into the ICM. As shown in dWhite et al. 2002t 
iBhattacharva et al. 20081: iBattaglia et al. 20101) . these effects 
change the SZ power spectrum by 10-30%. One possible 
way to include the uncertainty in gas ph ysics is the semi- 
analytic approach of modeling gas p hysics dBode et al. 2007: 
iBode et al. 20091: IShaw et al. 2010b and calibration the pa- 
rameters using X-ray observations of individual clusters. The 
other source of unc ertainty is the non -gaussian nature of the 
SZ power spectrum dShaw et al.2009|) . All these effects needs 
to be accounted for precision modeling of the power spec- 
trum. This will be essential to harness the full merit of the 
upcoming datasets. In the current st udy we choose to use the 
model by (|Komatsu & S eljak 2002|) to assess the impact of 
dark energy on the SZ power spectrum. 

4. RESULTS 

4.1. Cluster counts 

We consider two observational campaigns: upcoming X- 
ray survey eROSITA, and ongoing South Pole SZ survey. The 
first one, due in 2012, will cover almost half a sky (f^ w 
0.49), and have flux limit of F min = 3.3 x 10~ 14 ergs _1 cm~ 2 in 
energy band [0.5, 5.0] keV. South Pole Telescope will scan 
4000 square degrees (f s k y ~ 0.1) with limit of S m i n w 5mJy at 
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Fig. 4. — Expected redshift distribution of galaxy clusters for eROSITA 
Telescope. In black is our fiducial ACDM model, red and blue are EDE1, and 
EDE2, respectively, and in dashed lines we show results of EDE models when 
perturbations are neglected. Shaded regions are statistical Poisson errors plus 
10% systematic error coming from non-universality of the mass function. 
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Fig. 5. — The same as Fig. [4] but for the South Pole Telescope. 



frequency = 150Hz. Usin g mass-observable scaling rela- 
tions described in section l3Jl we can find minimum detectable 
mass for each survey, as well as redshift volume element be- 
ing observed (Figs. [2] and [3]). For the theoretical prediction 
of cluster abundance, we use Tinker et al. (2008) mass func- 
tion for the overdensity of 200p c (A « 778 in their notation), 
and we use their fitting formulas for the parameters of f(a) 
as the function of log A. Total number density of objects with 
masses above the given threshold as predicted by Tinker et 
al. (2008) is shown in lower left panels of Figures [2] and[3l for 
eROSITA and SPT survey, respectively. 

Figures 01 and \5\ show expected redshift distribution of 
galaxy clusters for the cosmologies considered here, where 
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we have taken redshift bins of Az = 0.1. As expected, SZ 
survey can detect meaningful number of clusters up to higher 
redshift, but has fewer objects at lower redshifts, due to the 
smaller sky coverage. Shaded area shows Poisson errorbar 
plus a systematic 10% uncertainty which describes our confi- 
dence in universal formula of the mass function. Dashed lines 
show expectations for EDE models if perturbations would be 
neglected. We see that presence of perturbations significantly 
affect growth of structure in the Universe, effectively reduc- 
ing erg . Since amount of this reduction is redshift dependent, 
it is not possible to mimic it with ACDM cosmology with dif- 
ferent erg (i.e. it is possible to mimic it at a particular redshift, 
but not at a redshift range). Most importantly, we see that per- 
turbations cannot be neglected when considering predictions 
from dynamical EDE models. 

Finally, Figures @] and [5] showcase that future large-scale 
structure results will be able to constrain EDE models signifi- 
cantly better than is currently the case. Both EDE1 and EDE2 
models, which represent models which are (a) indistinguish- 
able from ACDM at present (EDE1), or (b) at the upper limit 
of wo (EDE2), are distinguishable from ACDM cosmology 
with high significance (2a or more). Constraining power is 
larger on the EDE1 model which behaves the same as ACDM 
today but has a later transition from EDE like behavior (at 
redshift z — 5), as compared to the EDE2 model which is dif- 
ferent from ACDM today but has a much earlier transition 
from EDE like behavior (at redshift z = 9). Thus the cluster 
counts may constrain not only the current equation of state 
of dark energy, but also the redshift at which transition from 
EDE behavior occurs. 

4.2. SZ power spectrum 

In this section, we show the impact of the dark energy 
perturbation on the SZ power spectrum. The presence of 
perturbations in the dark energy changes the transfer func- 
tion, volume factor and the growth function of the universe, 
consequently the abundance of the halos changes. Since SZ 
power spectrum is proportional to the total abundance of the 
halos that have formed in the universe and volume of the 
universe, it depends strongly on the perturbations of dark 
energy. The impact of different early dark energy models on 
the SZ power spectrum is shown in figure [6l To assess the 
merit of the upcoming datasets to constrain early dark energy 
we assume two fiducial surveys: 

1) Planck like full sky survey: We assume 75% of the full 
sky coverage to exclude contaminations due to galactic fore- 
grounds, 7' resolution and 7 fiK per pixel as the sensitivity. 

2) ACT/SPT like higher resolution surveys: We assume 
a coverage of 4000 deg 2 , 1.2' angular resolution and a 
sensitivity of 20 fiK per pixel. 

The error on the SZ power spectrum for closely separated 
bins of size Al in / space can be written as 
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where f s k y is the fraction of the sky covered by the as- 
sumed surveys, w" 1 = [cr pix 6f whm /T C Mb] 2 , cr pix is the noise per 
pixel, 6f W hm is the resolution at full width at half maximum 
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Fig. 6. — The impact of dark energy perturbations on the SZ power spec- 
trum. Different early dark energy models are explained in section [2] The 
error bars represent a Planck like full sky survey and ACT/SPT like higher 
resolution survey. The error bars include cosmic variance including the non- 
gaussian contribution from the SZ trispectrum and the resolution of the sur- 
veys. We assume perfect removal of astrophysical contaminations from in- 
frared and radio galaxies and primary CMB. 

assuming a gaussian beam, T C mb is the CMB temperature 
( Jungman et al. 19961) . 
The trispectrum 7]/ (Poisson term) is given by 



Q = g 



dlnM **¥4xM,z,ef (17) 



dlnM 

As shown in fig. [H the upcoming measurements of SZ 
power spectrum will be able to distinguish between ACDM 
and different early dark energy models (that are not ruled out 
by current observations) with high significance. As in the case 
of the galaxy cluster counts, the SZ power spectrum is also 
sensitive to the redshift of transition from EDE-like to present 
day ACDM-like behavior. Current constraint s on the red - 
shift and width of transition are z t > 4, A, <, 0.2 (lAlam 20101) . 
while future observations of SZ power spectrum may even 
rule out EDE models with z t > 9, A t < 0.1 (EDE2). These 
results are of course for a fixed value of other non-dark en- 
ergy parameters, and we expect degeneracies with these in a 
full analysis. In the next section we attempt to constrain the 
EDE parameters using future simulated SZ power spectrum 
measurements and cluster counts. 

4.3. Likelihood analysis 

In this section we study the effect of the degeneracies be- 
tween the different cosmological parameters for EDE mod- 
els. We simula te the data according to WMAP7+BAO+# 
ACDM model dKomatsu et al 20 09). We generate the pri- 
mary CMB TT power spectrum simulated to match the Planck 
survey resolution and sky coverage. We also simulate the SZ 
power spectrum data as obtainable from ACT/SPT sur veys 
with the survey parameters discussed in detail in section l4~2l 
For studying the cluster data, we simulate SPT survey, assum- 
ing Tinker et al. (2008) number density of clusters; the details 
are given in section I47T1 The results for eRosita are expected 
to be similar. 

Typically, the cluster data alone is not enough to break 
the degeneracies between different parameters and we need 
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Fig. 7. — Confidence levels for a combination of primary CMB power spectrum from Planck survey + cluster counts obtainable from ACT/SPT survey. 





Fig. 8. — Confidence levels for a combination of primary CMB power spectrum from Planck survey +SZ power spectrum obtainable from ACT/SPT survey. 



the constraints from the primary CMB power spectrum. We 
add CMB data in the form of scalar Cl's simulated as per 
Planck specifications. We also apply a Gaussian prior on 
Ho as Ho = 70.4 ±3.6 km/s/Mpc, where the error is consis- 
tent with the recent results of the Hubble constant from the 
SHOE S (Supernovae an d //o for the Equation of State) pro- 
gram ( Ri ess et al. 2 009). We do not add any other possible 
future datasets such as the supernovae or Baryon Acoustic Os- 
cillations. 

In analyzing the data, we keep w m > -0.1 since we are in- 
terested in constraining models which have sufficient amount 
of dark energy at early times. The true ACDM model can- 
not be exactly reproduced by the EDE class of models since 
w m / -I, but in the limit wo = -l,a t —> (z t — » oo), A t — » 
these models replicate ACDM-like behavior. 

As shown in figure [71 CMB + cluster count anal- 
ysis can constrain the EDE parameters to wo< - 
0.9, a t < 0.2(z t > 4), A t < 0.3, which implies that cluster 
counts + CMB can constrain any deviation from ACDM be- 
havior upto at least z=4. For the SZ power spectrum data, 
we find that the EDE parameters are constrained to wo <, - 
0.85,0, < 0.15fc > 5.5), A t < 0.15 (see figure©. SZ power 
spectrum provide complementary information about dark en- 
ergy parameters compared to cluster counts. For example, 
SZ power spectrum put tighter constraints on the dark energy 
transition parameters compared to cluster counts. The matter 
density is reasonably reconstructed for both the cluster counts 
and SZ power spectrum. 

5. CONCLUSIONS 



In this work, we have studied early dark energy models in 
light of future galaxy cluster data. We study two EDE mod- 
els, both of which are allowed within the current observational 
constraints. The first model (EDE1) transits from EDE behav- 
ior at redshift z > 5 to ACDM-like behavior at present, while 
the second model (EDE2) transits at a higher redshift (z — 9) 
to an equation of state w>o = -0.9, which is close to but not 
identical to ACDM. For both the models, as well as for the 
fiducial ACDM cosmology, we show the predictions for clus- 
ter abundance and SZ power spectrum. For that purpose, we 
consider two instruments - eROSITA for X-ray surveys, and 
ACT/SPT for SZ effect. We also do a likelihood analysis of 
data simulated replicating the SPT and eRosita survey specifi- 
cations, along with the simulated Planck CMB data, to obtain 
the constraints on the EDE parameters from future surveys, 
and find that future galaxy cluster counts and SZ power spec- 
trum can put competitive constraints on these parameters. 

It is worth noting some apparent differences between our 
work here, and some other recent works on EDE signatures. 
For example, Fedeli et al. (2009) also considered several EDE 
models, but their cluster counts are higher than for the fidu- 
cial ACDM model, although ag for EDE models are com- 
monly lower in their work (as it is the case here as well). That 
is because for their choice of EDE models Hubble parame- 
ters, E(z) differs significantly from the corresponding ACDM 
model. This affects minimum observable mass as a function 
of redshift; their EDE models have much lower mass thresh- 
old than ACDM, thus resulting in higher observable cluster 
counts. In contrast, we have chosen models which have E(z) 
similar to that of the corresponding ACDM model, to high- 
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light the difference coming from the perturbations rather than 
the background expansion. Similar difference from this work, 
due to the radical departure of E(z) from ACDM, can be seen 
in Waizmann & Bartelmann (2009), for the effects of EDE 
on SZ power spectrum. In some other cases differences can 
arise depending whether one uses low-k, CMB normalization 
of the power spectrum as done here, or ag normalization as 
in Sadeh et al. (2007), although it is clear that a successful 
model has to fulfill both constraints. 

We show that the inclusion of dark energy perturbations 
has a major effect on the matter power spectrum, therefore 
increasing the possibility of discriminating EDE models from 
ACDM using large scale structure probes. Neglecting the dark 
energy perturbations leads to severe underestimation of the 
imprint which the sharp transition in dark energy equation of 
state leaves on the dark energy perturbations and therefore on 
the matter power spectrum. We show that the models con- 
sidered here which are allowed by the current observations 
can be ruled out using the future galaxy cluster probes. It is 
also interesting that both the cluster counts and the SZ power 
spectrum are sensitive to the redshift at which the transition 
between early and present day dark energy occurs. We expect 
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to put strong constraints on the equation of state of dark en- 
ergy at present using low redshift geometric observables (such 
as the luminosity distance of type la SNe and Baryon Acous- 
tic Oscillations peaks). These probes of geometry of the Uni- 
verse, however, are insensitive to the high redshift behavior of 
dark energy. The galaxy cluster observables such as their red- 
shift abundance and the SZ power spectrum, although at low 
redshift, are sensitive to the perturbations in the early dark en- 
ergy models, and hence would be able to put constraints on 
the redshift of transition from EDE behavior to present-day, 
ACDM-like behavior. With the ongoing and future cluster 
surveys in microwave, optical and X-ray waveband, galaxy 
clusters will be able to provide strong constraints on the dy- 
namical dark energy sector. 
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